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ABSTRACT 

Highly unequal-mass ratio binaries are rare among field brown dwarfs, with the mass 
ratio distribution of the known census described by qi( 4 - 9=1=0 - 7 ) . However, such sys- 
tems enable a unique test of the joint accuracy of evolutionary and atmospheric mod- 
els, under the constraint of coevality for the individual components (the "isochrone 
test"). We carry out this test using two of the most extreme field substellar bi- 
naries currently known, the T1+T6 e Ind Bab binary and a newly discovered 0.14" 
T2.0+T7.5 binary, 2MASS J12095613-1004008AB, identified with Keck laser guide 
star adaptive optics. The latter is the most extreme tight binary resolved to date 
(q ~ 0.5). Based on the locations of the binary components on the H-R diagram, 
current models successfully indicate that these two systems are coeval, with internal 
age differences of \og{age) = -0.8 ± 1.3 (— 1.0±JJ) dex and 0.5±£| (0.3±g;|) dex for 
2MASS J1209— 1004AB and e Ind Bab, respectively, as inferred from the Lyon (Tuc- 
son) models. However, the total mass of e Ind Bab derived from the H-R diagram 
(~80 Mj up using the Lyon models) is strongly discrepant with the reported dynamical 
mass. This problem, which is independent of the assumed age of the e Ind Bab system, 
can be explained by a «50-100 K systematic error in the model atmosphere fitting, 
indicating slightly warmer temperatures for both components; bringing the mass de- 
terminations from the H-R diagram and the visual orbit into consistency leads to an 
inferred age of ~6 Gyr for e Ind Bab, older than previously assumed. Overall, the two 
T dwarf binaries studied here, along with recent results from T dwarfs in age and mass 
benchmark systems, yield evidence for small (~100 K) errors in the evolutionary mod- 
els and/or model atmospheres, but not significantly larger. Future parallax, resolved 
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spectroscopy, and dynamical mass measurements for 2MASS J1209— 1004AB will en- 
able a more stringent application of the isochrone test. Finally, the binary nature of 
this object reduces its utility as the primary T3 near-IR spectral typing standard; we 
suggest SDSS J1206+2813 as a replacement. 

Subject headings: binaries: general, close — stars: brown dwarfs — infrared: stars - 
techniques: high angular resolution 



1. Introduction 



Wide-field surveys have now identified nearly a thousand ultracool dwarfs (spectral types 
later than M6) in the solar neighborhood, spanning temperatures and (inferred) masses that range 
from the stellar/substellar boundary defined by the hydrogen-burning minimum mass (~75 Mj up ), 
down to the proposed substellar /planetary boundary demarcated by the deuterium-burning limit 
(«13 Mj up ). This census has enabled extensive characterization of the spectrophotometric prop- 
erties of these low-luminosity, low-temperature objects. However, studies of the substellar field 
population are inevitably hindered by the unknown properties of individual objects, making it 
difficult to disentangle the effects of varying ages, masses, and compositions on the underlying 
physics. In principle, a sample of star clusters with complete membership down to very low masses 
and encompassing a range of ages would be ide al. However, such a sensitive census onl y exists 
for a handful of young (~few Myr) clusters (e.g. lLodieu et al.l 120071 : iLuhman et al.ll2007al ). where 
the lowest mass members are easiest to detect. With increasing age, dynamical evolution leads 
to depl etion of the lowest mass members such that older clusters may be largely devoid of brown 
dwarfs (jAdams et alJbood ; iBouvier et alJbood ). 



One avenue for circumventing these limitations is through the study of binary brown dwarfs, 
as binaries constitute systems of common (albeit unknown) age and metallicity. About 15% of 
field brown dwarfs are tight binary systems, resolved largely by Hubble Space Telescope or ground- 
based adaptive optics (AO) i maging and composed mostly of nearly equal-luminosity (and thus 
equal-mass) components (e.g. iBurgasser et al.1 120071 ) . These "mini-clusters" have proven to be 
fertile grou nd for understanding the physical processes that govern the emergent spectra of brown 
dwarfs (e.g. lLiu fc Leggettll2005l : LLiu et al.ll2006l : IBurgasser et al.ll2006dr) and for testing theoretical 



models with dynamical mass determinations (e.g 



Liu et al.l 1200a : iDupuv et alj|2009d ; iKonopackv et al 



atero Osorio et al 



Zapatero 
alJ boid : 



Dupuv et al 



2004 



201(1 ) . 



Ireland et al J 1200a 



As part of our ongoing program to study the multiplicity and physical properties of substellar 
binaries using laser guide star adaptive optics (LGS AO), we present here the discovery of the 
binarity of 2MASS J12 095613- 1004008. herein after 2MASS J1209-1004. This object was identi- 



fied in 2MASS data by IBurgasser et al.l (12004! ) and spectrally typed by them and by 



(|2006l ) as a T3 dwarf based on integrated-light near-IR spectroscopy. IBurgasser et al 



Chiu et al 



2006b|) de 



fine 2MASS J1209— 1004 as the primary T3 spectral type standard for their near-IR classification 
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scheme for T dwarfs. This object has not previously been targeted with high angular resolution 
imaging. What distinguishes 2MASS J1209— 1004AB from most previously known substellar field 
binaries is the large IR brightness difference between its two components, indicating an atypi- 
cal mass ratio compared to the plethora of nearly equal-mass binaries. As we describe below, 
2MASS J1209— 1004AB provides a new opportunity to test theoretical models of brown dwarfs - 
by requiring the model-derived ages of the two components to be consistent with coevality of the 
system, an approach which we refer to here as the "isochrone test." 



2. Observations 

We imaged 2MASS J1209-1004 on 22 April 2007 and 16 January 20 08 UT using the sodium 



LGS AO system of th e 10-meter Keck II Telescope on Mauna Kea, Hawaii (jWizinowich et al. 



2006 



van Dam et ajj|2006j ). We used the facility IR camera NIRC2 with its narrow field-of-view camera, 
which produces a 10.2" x 10.2" field of view. Conditions were photometric for both runs. The 
LGS provided the wavefront reference source for AO correction, with the tip-til t motion measured 



contemporaneously from the R = 14.7 mag field star USNO-B1.0 0799-0230529 (IMonet et al.ll2003 ) 
located 68" away from 2MASS J1209— 1004. The LGS brightness, as measured by the flux incident 
on the AO wavefront sensor, was equivalent to a V ~ 9.6 — 10.0 mag star. In April 2007, we 
obtained images with the MKO J (1.25 /xm), H (1.64 /xm), and K (2.20 /jm) filters. In January 
2008, we obtained images with the CH^s filter, which has a central wavelength of 1.592 //m and a 
width of 0.126 /im. 



The images were analyzed in the same fashion as our previous LGS papers (ILiu et al.l 12008 



20061 ). The raw images were reduced using standard methods of flat-fielding and sky-subtraction. 
The binary's flux ratios and relative astrometry were derived by fitting an analytic model of the 
point spread function (PSF) as the sum of three elliptical gaussians. The measurements were 
made on the individual images, with outlier images of much poorer quality excluded. The averages 
of the fitting results were adopted as the final measurements and the standard deviations as the 
errors. The final flux ratios were used to correct the Strehl ratio measurements of the images 
for the contamination from the light of the secondary component. We did not correct the relative 
astrometry for instrumental optical distortion, as the size of the effect as predicted from a distortion 
solution by B. Cameron (priv. comm.) is insignificant compared to our measurements uncertainties. 

In order to validate our measurements, we created myriad artificial binary stars from LGS 
images of single stars with comparable Strehl and FWHM as the science data. For each fil- 
ter, our fitting code was applied to artificial binaries with similar separations and flux ratios as 
2MASS J1209— 1004AB over a range of PAs. The simulations showed that the random errors are 
reasonable and any systematic offsets are small. In cases where the measurement errors from 
the artificial binaries were larger than those from the 2MASS J1209— 1004AB measurements, we 
conservatively adopted the larger errors. 
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We adopted a pixel scale of 9.963±0.005 mas /pixel and an or ientation for the detector's +y axis 
of +0.13 ± 0.07° for NIRC2 as determined bv lGhez et all (j2008h . We computed the expected shift 
in the relative astrometry of the two component s due to differential chromatic refraction (DCR) 
in the same manner as in our previous work (e.g. iDupuv et al.ll2009d ). The DCR effects are small 
(~l-3 mas) compared to the measurement errors, and therefore we did not account for them. 
Table Q] presents our final Keck LGS measurements, and Figure Q] shows our Keck LGS images. 



3. Results for 2MASS J1209-1004AB 

3.1. Proof of Companionship 

The near-IR broadband colors of 2MASS J1209— 1004B are extremely blue, which provides 
strong circumstantial evidence that it is a p hysically associated late-T dwarf companion and not 



a background object. iBurgasser et al,l (|2004l ) measure a proper motion for 2MASS J1209— 1004 of 
0.46 ± 0.10"/yr at a PA of 140 ± 8°. Considering this proper motion, examination of the Digitized 
Sky Survey images shows no plausible optical counterpart if component B were a background 
object. Our two epochs of Keck LGS imaging separated by 9 months are consistent with no change 
in the relative position of the two components. If the companion were a background object, in 
Jan 2008 the system should have had a separation of 0.48 ± 0.07" at a PA of 317 ± 6°, which 
is 4.8(7 discrepant with the observed separation, neglecting the effect of parallax. Finally, our 
LGS photometry indicates the presence of .ff-band methane absorption in component B (§ I3.2[) . 
inconsistent with the background star hypothesis. Thus both the astrometry and photometry 
independently indicate that the system is a true physical binary. 



3.2. Resolved Photometry 



We use our measured flux ratios and the published JHK photometry from IChiu et al.l (120061 ) 
to derive resolved IR colors and magnitudes for 2MASS J1209— 1004AB on the MKO system. For 
the CH^s data, we synthesize i n tegrat ed-light photometry of 15.05 ± 0.03 mag from the near- 
IR spectrum of IBurgasser et al.l (120041 ). flux-calibrated to the broadband MKO photometry of 
H = 15.24 ±0.03 mag. We assume that the components of 2MASS J1209-1004AB are themselves 
single, not unresolved binaries. Table [2] presents the resolved photometry. Comparing to known 
ultracool dwarfs, Figure [2] shows that component A has IR colors most typical of T2-T3 dwarfs. 
The uncertainties are larger for component B, but its near-IR colors appear to be comparable or 
even bluer than the very latest T dwarfs. 

The resolved (CH4S — H) colors provide more accurate estimates of the spectral types than the 
broadband colors, as t he ff-band methan e absorption correlates well with overall near-IR spectral 
type (e.g., Figure 2 of lTinney et al.ll2005l ). We convert these colors to near-IR spectral type using 
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the polynomial fit from lLiu et al.1 (|2008l ): 

SpT = 19.40 - 19.698 x (CH 4 s - H) - 8.3600 x (CH iS - H) 2 



where SpT = 20 for TO, = 21 for Tl, etc. The RMS scatter about the fits is 0.3 subclasses. The 
fit was originally derived for objects from TO to T8, but as discussed below the more recently 
discovered T9 objects are still well-described by this relation. 

The observed {CH 4 s - H) colors for 2MASS J1209-1004AB give spectral types of T1.6 ± 0.9 
and T9 ± 2 for components A and B, respectively, where the spectral type uncertainties come 
from propagation of the errors in the colors. Note that the (CH^s — H) color of component B 
(—0.8 ± 0.3 mag) is bluer than any of the objects used to define the polynomial fit, so the spectral 
type estimate represents a small extrapolation of the relation. For comparison, the colors of the 
T8.5 dwarfs ULAS 1238+0953 and ULAS 2146-0010 are -0.65 and -0.63 mags, respectively, 
and those of the T9 dwarfs ULAS 0034-0052, CFBDS 0059-0114, and ULAS 1335+1130 are 
0.56, —0.63, and —0-66 mags, respectively, as synthesized from their published near-IR spectra 



([Warren et al.ll2007l ; iDelorme et al.ll2008l ; iBurningham et al.ll200a 12009 ) . (These very late-T dwarfs 
are not included in the Liu et al. fit, but their inclusion would make no significant difference in the 
polynomial relation.) Overall, component B shows ff-band methane absorption comparable to the 
latest known T dwarfs. 



3.3. Composite Spectra and Resolved Spectral Types 



We modeled the integrated-light spectrum of 2MASS J 1209— 1004 as the composite sum of 
two template T dwarfs, as an additional means to determine the resolved spectral types. The 
templates were chosen from available near-IR spectra for dwarfs that have similar colo rs to those of 



2MA SS J1209-1004 A and B. We produced variou s fits using SDSS J0758+324 7 fT2:lKnapp et al 



(T8; IBurgasser et al. 



(T7.5; IPinfield et al. 



20041 ). SDSS J1254-0122 (T2; |Leggett et al]l200d ) an d SDSS J1521+0131 ( T2; Knapp et alJl2004 ) 
as te mplates for the compo nent A and Gl 570D fT7.5: feurgasser et al.lhoooh'l . 2MASS J0415-0935 



2002h . ULAS J1017+0118 CT8p: IBurningham et alJ200a ) and ULAS J1315+0826 



200a ) as templates for component B. The JHK colors of the template dwarfs 



and the resolved binary are within 0.2 mag for the component A and 0.4 mag for the component B 
(whose observed photometry is much more uncertain). 

Before combining each template pair, we scaled each template spectrum to a distance of 10 pc. 
Trigonometric parallaxes are known for SDSS J1254-0122 and 2MASS J0415-0935. For the other 
dwarfs dis tances were estim ated using the relationships between absolute magnitudes and spectral 
type from lLiu et al.1 (|2006l ). We explored fits using both the "bright" and "faint" relations, which 
differ significantly for early-type T dwarfs. For each pair, after shifting the flux to lOpc, the flux 
at each wavelength was summed, and the result compared to the observed unresolved spectrum 
for 2MASS J1209-1004AB. There is no parallax measurement for 2MASS J1209-1004AB, and 
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so we shifted the observed spectrum to that of the synthesized spectrum, determined the implied 
distance, and compared the resulting spectra. 

In the end, model composite spectra composed of a T2 primary and T7.5 secondary showed 
good agreement with the observed integrated-light spectrum (Figure [3]). Any combination using the 
peculiar T8 dwarf ULAS J1017+0118 failed to reproduce the observed strength of the CH4 feature 
at 1.6 fim. The best fits were obtained using the T2 dwarf SDSS J1254— 0122 in combination with 
either of the T7.5 dwarfs Gl 570D or ULAS J1315+0826, or using the T2 dwarf SDSS J1521+0131 
in combination with the T8 dwarf 2MASS J0415— 0935. For the last combination to work, the 
distance for SDSS J1521+0131 must be derived using the "faint" Liu et al. photometric relation. 
The implied distance for 2MASS J1209-1004AB is 23 pc from the SDSS J1254-0122 fits, or 18 pc 
from the SDSS J1521+0131 fit. 

Combining the results from the JHK and {CH^s — H) colors and this spectral modeling, we 
adopt a spectral type of T2 ± 0.5 for 2MASS J1209— 1004A, as all the available estimates are in 
good agreement. 

The typing for component B is more challenging, given its larger photometric errors and its 
relative faintness (and thus modest contribution to the integrated-light spectrum). In the spectral 
modeling, T7.5 templates provided the best matches. Our attempts to use T8 dwarfs failed to pro- 
duce a good match, though this was possibly limited by the small number of templates with compa- 
rably blue near-IR colors as 2MASS J1209-1004B (2MASS J0415-0935 and ULAS J1017+0 118). 



The polynomial fits for absolute magnitude as a function of spectral type from lLiu et al.l (|200Q ) give 
M{J) = {15.1, 15.5, 15.9, 16.3} mag, M(H) = {15.5, 15.9, 16.3, 16.7} mag, and M(K) = {15.6, 
16.0, 16.4, 16.7} mag for near-IR spectral types of T6.5, T7, T7.5, and T8, respectively. Based 
on the photometric distance to component A (§ \3.5\i . the absolute magnitudes of component B are 
M( J, H, K) = 15.7 ± 0.4, 16.5 ± 0.5, 16.9 ± 0.7 mag, in good agreement with T7.5. Given the con- 
cordance between this and the spectral modeling, we adopt a spectral type estimate of T7.5 ± 0.5 
for component B. 



3.4. T3 Near-IR Spectral Type Standard 



Burgasser et al.l (|2006bl ) chose 2MASS J1209— 1004 as the primary T3 spectral standard for 



near-IR classification of T dwarfs. Its binary nature now lessens its utility for this purpose, espe- 
cially at the bluer wavelengths where the contribution of the secondary to the integrated-light spec- 
trum is greater. The secondary standard proposed by Burgasser et al. was SP SS J1021— 0304AB, 
which also has been resolved into a binary (T1+T5; iBurgasser et al.1 l2006dl ) . For both sources, 
their binary nature make them unappealing as the T3 spectral standard, since the secondary 
contributes significantly to the integrated near-infrared light (e.g., Figure 3)0 Furthermore the 



The mid-infrared regime becomes increasingly dominant for cooler T dwarfs (e.g. Figure 3 of Cushing et al. 
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primary of neither system is of T3 type. Options for a different spectral standard are limited, due 
to the modest number of T3 dwarfs known. At the relatively brighter magnitudes, the choices are 
SDSS J120602.51+281328.7 (T3; J = 16.5 mag) and SPSS J141530 .05+572428.7 (T3±l; J = 16.7), 
both of which were identified and classified by IChiu et al.l (|2006l ). SDSS J1206+2813 appears to 
be single in 0.06" FWHM images obtained by us with Keck LGS (Liu et al., in prep.) while no 
high angular resolution imaging is available yet for SDSS J1415+5724. As the individual spectral 
indices for SDSS J1206+2813 are more consistent with a T3 typing and it appears to be single, we 
suggest that it be used as the T3 spectral standard. 



3.5. Photometric Distance 



With an estimated spectral type of T2.0 ± 0.5 for component A, we use the iLiu et al.l (|2006l ) 
polynomial fits of absolute magnitude as a function of spectral type to determine a photometric 
distance. We use the J-band photometry; the H and A'-band data produce nearly identical results 
but with somewhat larger uncertainties. Errors in the photometric distance include the uncertain- 
ties in the apparent magnitude (0.05 mag), the uncertainty in J-band absolute magnitude arising 
from the spectral type uncertainty (0.05 mag), and the intrinsic dispersion about the fit (0.4 mag). 
Liu et al. provide both a "bright" and a "faint" relation, depending on the inclusion or exclusion of 
candidate overluminous objects of early/mid-T spectral type, respectively. The two relations give 
consistent answers within the errors, 24 ± 4 and 18 ± 4 pc, respectively. We average to arrive at a 
final photometric distance of 21 ± 4 pc0 

Since the spectral templates used for modeling the composite spectra have distance determi- 
nations (§ 13. 3p . the scaling factor used to match to 2MASS J1209— 1004AB also provides a distance 
estimate to the binary, with a value of 18 or 23 pc f§ 13.31). An e rror estimate from this approach is 
uncertain. The T2 dwarf SDSS J1254— 0122 (jLeggett et al.ll200C ) provides a well- matching template 
for component A and has parallax with an uncertainty of 4% (IVrba et al.ll2004l ). which represents 



the ve ry smallest possible uncertainty. At the other extreme, the photometric relations of lLiu et al 



(|2006l ) have a RMS scatter about the fit of 15-20% in distance. Adopting a distance error of 10% is 
a reasonable compromise; for instance, Liu et al. demonstrate that fitting the composite spectrum 
of the T1+T6 binary e Ind Bab produces a distance estimate that agrees to 10% with the measured 
parallax. Thus, the distance estimates from the spectral modeling agree with the photometric 
distance of 21 ± 4 pc. We conservatively adopt the latter in our subsequent analysis^ 



2006) and so the T7.5 2MASS J1209— 1004B is also expected to make a significant contribution to the integrated 
2MASS J1209-1004 flux at ~10 /itm. 

2 In a similar fashion, we find a photometric distance to component B of 19 ± 5 pc based on its type of T7.5 ± 0.5. 
This value is consistent with the distance derived component A, but it cannot be considered truly an independent 
measurement, since we use the photometric distance in estimating the spectral type for component B in § 13.31 

3 The slightly larger distance from the spectral modeling results from the fact that SDSS J1254— 0122 is slightly 
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3.6. Bolometric Luminosities 

3.6.1. Combined System 
There is no published value for the integrated-light bolometric luminosity of 2MASS J1209— 1004, 



so we computed it dir ectly from its near-IR spectrum (jBurgasser et al.ll2004j ) . flux calibrated with 



the IChiu et al.l (|2006l ) photometry, and combined with mid-IR photometry. For the latter, IRAC 
photometry is available in the Spitzer archive, obtained as part of the Cycle 1 GTO program 35 on 
UT 2005 June 13. Exposure times were 30 s, and a 5-position small dither pattern was used. The 
binary is unresolved in the IRAC images, and an 8" aperture was used for the photometry. We 
determined magnitudes of [3.6] = 14.02 ± 0.004, [4.5] = 13.49 ± 0.004, [5.8] = 13.33 ± 0.016, and 
[8.0] = 13.06 ± 0.06 mags. A 0.03 mag error should be added in quadrature to the quoted random 
errors to account for systematic effects. 

To derive the integrated-light Lboi, we numerically integrated the near-IR spectrum and the 
mid-IR photometry, interpolating between the gaps in the data, extrapolating at shorter wave- 
lengths to zero flux at zero wavelength, and extrapolating beyond 4.5 (im assuming a blackbody. 
We determined the luminosity error in a Monte Carlo fashion, by adding randomly drawn noise to 
our data over many trials and computing the rms of the resulting luminosities. In this process, we 
accounted for the noise in the spectrum, the errors in the MKO photometry used to flux-calibrate 
it, and the errors in our mid-IR photometry estimates. Before accounting for the error in the photo- 
metric distance, we found a total bolometric luminosity of log(Lb i/L©) = —4.61 ±0.02 dex. With 
the 20% uncertainty in the photometric distance included, the Lboi uncertainty becomes 0.15 dex. 



3.6.2. Individual Components 



To compute the bolometric luminosities (Lboi) for each component, we apply bolometric cor- 
rections to the resolved near-IR photometry based on the compilation in Appendix [A] Despite the 
relatively larger RMS about the BC fit at J-band, the photometry of the component B is most 
accurately determined in this bandpass, so we use it to to determine Lboi and associated quantities. 
The results are given in Table [2j The Lboi uncertainties are composed of the uncertainties from the 
apparent magnitudes, from the bolometric corrections arising from the spectral type uncertainty, 
the intrinsic dispersion about our polynomial fits, and the photometric distance. The last item 
dominates the uncertainties in the Lboi values. (However, the luminosity ratio of two components 
is known to much higher precision than the individual Lboi values, since the ratio is independent 



brighter than the values from the adopted Liu et ah polynom ial fi ts. In fact, it has been suggested this system may 
be a binary based on this modest over luminosity |Liu et alJbood : lBurgaase3l2007h. though it h as not been resolved 
as such in high angular resolution imaging from HST or Keck LGS AO jBurgasser et al.ll2006dl ; Liu et al. in prep.) 
and modeling of its integrated-light s pectra suggests the sy stem is more likely be a relatively young, single object as 
opposed to a near-equal mass binary ( Stephens et al.ll2009T ). 
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of the distance uncertainty.) As a consistency check, we note that the sum of the individual com- 
ponent luminosities derived from BCj (—4.51 ±0.18 dex) agrees well with the direct computation 
from the integrated light spectrum (—4.61 ± 0.15 dex). 



3.7. Physical Properties from Evolutionary Models 



With the Lbol determinations for the individual components and an assumed ag e, evolutionary 



mode ls can then fully determine the physical propertie s of the two componen ts (e.g. ISaumon et al 



2000J). Figure H] shows the results graphically using the lBurrows et al.l (j 19971 ) evolutionary models, 
where we have adopted a generic age range of 0.1-10 Gyr for the 2MASS J1209— 1004AB system. 
The figure indicates the complete possible set of {T e fj, log(<?)} values for the two components, based 
on our Lbol determinations. 

Table [3] summarizes the physical properties derived from this approach, for a range of assumed 
ages. Note that for computing relative quantities between the components (e.g., luminosity ratio 
and mass ratio), we take care to account for the covariance of the measurements due to the fact 
that the uncertainty in the distance is common to both components, and thus uncertainties in 
these relative quantities are smaller than for the absolute value of these quantities for the individual 
components. To do this, we draw values for all the relevant measurements in a Monte Carlo fashion, 
keeping track of all the calculations of the physical parameters in every trial, and then compute 
the final result and RMS from the ensemble of trials. 



To compute the estimated orbital period of the system, following iTorresI (Il999l ). we assume 
random viewing angles and a uniform eccentricity distribution from e = — 1. This gives a 
multiplicative correction factor of LlO^gg (68.3% confidence limits) for converting the projected 
separation into a semi-major axis. The estimated orbital periods range from 16-26 years for ages of 
0.5, 1.0, and 5.0 Gyr, albeit with significant uncertainties. Dynamical masses for ultracool binaries 
are possible fr om orbital monitoring covering >30% of the orbital period (e.g. iBouv et al.l 12004 ; 
Liu et al.ll2008l ). Thus, 2MASS J1209— 1004AB warrants continued high angular resolution imaging 
to determine its visual orbit, combined with a parallax measurement to compute its total mass. 



4. Discussion: The Isochrone Test 

The 2MASS J1209— 1004AB system is a rare example of an ultracool binary with a highly 
unequal mass ratio (q ~ 0.5), as most ultracool binaries are composed of nearly equal components 
(Table . This novelty is highlighted in Figure which shows the estimated mass ratios of all 
known field brown dwarf binaries, where we choose field systems in which the pr imary has an esti 



mate d spectral type of L4 or later as a proxy for the stellar /substellar boundary (jKirkpatrick et al 



19991 ) . As explained below, we specifically focus on ultracool binaries composed of two brown 



dwarfs, as opposed to those that contain one or two very low-mass stars. 
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Under the conservative assumption of coevality, we can use 2MASS J1209— 1004AB to test 
current theoretical models of brown dwarfs, by assessing whether the two components lie along a 
single isochrone on the Hertzsprung-Russell (H-R) diagram. We essentially treat the system as a 
star cluster composed of two members, with an unknown age but requiring that it is the same for 
components A and B. Hereinafter, we refer to this approach as the "isochrone test." 

Placement on the H-R diagram requires Lbol an d T e Q. The former can be well-determined for 
ultracool dwarfs (Appendix lAl and references ther ein). However, robu st temperature determina- 
tions are still an active area of investigation (e.g. ICushing et al.ll2008l ). Atmospheric models for 
L and T dwarfs have not been confirmed by direct measurement of brown dwarf radii, only checked 
for consistency against evolutionary models using single brown dwarfs that are companions to 



stars of known age (e.g. ISaumon et al. 



2006 



Leggett et al.l 120081) and the few brow n dwarf bina- 



Liu et al. 


2008|: 


Dupuv et al. 


2009b 



ries with dynamical mass determinations (e.g. 
low-temperature atmospheres is challenging. The molecular opacity line list for CH4, a species 
important for both the L and T dwarfs, is incomplete or non-existent below 1.6 fim. The FeH line 
list is incomplete in the 1.2-1.3 fim region, which is important for L dwarfs. NH3 is important for 
very late-type T dwarfs, but there are no opacities calculated below 1.4 ftm. Also, calculation of the 
very strongly pressure- broadened K I resonance doublet (0.78 /xm) is difficult, and hence models 
do not currently reproduce the far-red region of T dwarf spectra very well. Finally, the treatment 
of the grains condensates is difficult, which impacts the 1.0-1.6 [im region of ~L3 to ~T2 types. 

Thus, our current ability to carry out the isochrone test represents a test of the joint accuracy 
of evolutionary and atmospheric models, not of either specific class of modelslfl The fact that the 
isochrone test relies on components of very different masses also may make it a useful signpost for 
mass-dependent problems, which would not be revealed when studying the (much more common) 
nearly-equal mass binaries. Of course, the coevality requirement is widely applied in studies of 
star clusters. Its particular characteristic in the case of substellar binaries arises from the fact that 
brown dwarfs steadily cool over their lifetime, and thus the isochrones for substellar objects are 
more dispersed on the H-R diagram compared to stellar ,sochrones@ 



4 Of course, for this purpos e we cannot use temperatures that are derived from the evolutionary models themselves 
fe.g. lGolimowski et a l. 2004a), where the observed Lbol, an assumed age, and the evolutionary-model-derived radii are 
used to compute T e g. If we did, then any information gleaned from the H-R diagram positions would automatically 
agree with the evolutionary models, i.e., this would partly be circular reasoning. 

5 One limitation on the isochrone test is that the currently available evolutionary models do not t ypica lly incorpo- 
rate the same set of model atmospheres used in fitting observed spectra. However, IChabrier et al.l (|2000l ) show that 



the predicted luminosity evolution is mostly robust against the adopted model atmosphere. 

6 The isochrone test described here is conceptually identical to some previous studies of young star-forming regions. 
Although the underlying physical processes are different for pre-main sequence stars and field brown dwarfs, the 
luminosities of both classes of objects are strongly dependent on age. For young regions, coevality has been used to 
chec k the accuracy of evolutionary models and the spectral type-to-effective temperature scal e, using multiple systems 
(e.g. Hartigan et al.ll 1994 ; White et al. 19991 ; Schaefer et al. 20081 ; Kraus fc Hillenbrandlboolj ) and entire clusters (e.g. 



Luhman 



1999m . Note for such young objects, a complicating factor is the possible impact of accretion on the H-R 
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4.1. 2MASS J1209-1004AB 



We now subject 2MASS J1209— 1004AB to the isochrone test. A parallax measurement is not 
available, so we use the photometric distance estimate (§ 13. 5 p — this is not a limiting factor in the 
analysis. For the effective temperatures of each component, we consider previous analysis of the 
spectra for similarly typed T dwarfs. 

Three T2 dwar fs have been analyzed in such a f ashion: SPSS J0758+3247 (1100-1200 K ; 
Stephens et ailhood ). SPSS J1254-0122 ( 1100-1200 K: lstephens et al.ll2009l : Tcushing et alibood l. 



and HN Peg B (1115 Ki lLeggett et al.ll2008l ) II The model atmospheres were spaced by 100 K (except 
for HN Peg B), and the estimated systematic uncertainties are ~100 K due to, e.g., flux calibration, 
choice of wavelength region for fitting, and weighting scheme. Thus we adopt 1150 ± 100 K as the 
temperature for 2MASS J1209— 1004A. No Tl or T3 dwarfs have been s ubjected to such analy sis, 



but the T e g fitted for early/mid-T dwarfs are all relatively constant (e.g. IStephens et al.ll2009l ). in 



accord with the T e g measurements derived from Lhni data that the L /T transition for field objects 



occurs at nearly constant temperature (e.g. iGolimowski et al.ll2004al ). 



Four T7.5 dwarfs have been analyzed to date: HP 3651B (820-830 K; iLeggett et al 



2MASS J1114-2618 (725-775 K: lLeggett et all 120071). G l 570 P (800-821 K: ISaumon et al. 
and 2MASS J1217-0311 (850-950 LK- Tsaumon et al 



2007) 



2006) 



20071 ). To encompass the potential uncertainty 



in the spectral ty pe of component B, we also consider the T 8 dwarf 2MASS J0415— 0935 (725- 
775 K, 710 ±40 K 'lSaumon et al.ll2007i : burgasser et al.ll2008bh Fl No T7 dwarfs have independently 
determined T e g 'so Taking the unweighted mean and RMS of the results gives T e g- = 800 ± 70 K as 
our adopted value for 2MASS J1209-1004B. 

Figure [6] shows 2MASS J1209— 1004AB on the H-R diagram. We compute t he individual 



comp onent ages bas ed on their H-R diag ram positions relative to the Lyon/CONP (jBaraffe et al 



20031 ) and Tucson (IBurrows et al.1 119971 ) evolutionary models, using a Monte Carlo method to 
account for the uncertainties. For each component, we draw trial values for Lbol and T c g from 
a normal distribution, which are then used with finely interpolated tabulations of the models to 
derive a set of ages. We then take the difference of the component ages and compile the distribution 



diagram positions ( Baraffe et al. 20091 ) . an effect which is not relevant for old (field) brown dwarfs. 



7 The first two objects were analyzed with direct fitting of model atmospheres to their observed spectra. For 
HN Peg B, the model atmosphere matching the T c ff derived from evolutionary models was shown to be consistent 
with the observed flux-calibrated spectrum. 

8 For Gl 570D and 2M ASS J0415— 0935, the model atmosphere matching the T e g derived from evolutionary mod- 
els was shown to be consistent with the observed flux-calibrated spectrum, as opposed to directly fitting model 
atmospheres to the spectra. 



iBurgasser et al.l ( 2006aT ) have determined T e g's for a sample of mid/late-T dwarfs based on multi-band flux ratios 
measured from low-resolution near-IR spectra. However, these are calibrated against the T c g and log(g) of Gl 570D 
derived directly from evolutionary models. And thus these Tog's cannot be used to place objects on the H-R diagram, 
since this would be circular reasoning. 
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resulting from 10 6 trialscj Figure [7] plots the result, which is well-described by a gaussian. Table [5] 
gives the derived ages and age difference. At the la level, the models successfully indicate that the 
two components are coeval, with a log(a<?e) difference of — 0.8^'| and — l.O^}^ dex inferred from 
the Lyon and Tucson models, respectively. 



4.2. e Ind Bab 



We also apply the isochrone test to the e Ind Bab system, w hich is a wide (1459 AU) separation 
binary companion to the K5 V star e IndA fscholz et alJbood) The b inary is composed of a Tl 



and a T6 component (jMcCaughrean et al. 



2004 



Burgasser et al.ll2006bl ). with a -fT-band flux ratio 



of 2.18 ±0.03 mag. The system is quite similar to 2MASS J1209-1004AB, but with more precisely 
determined physical parameters and therefore offers a more stringent application of the isochrone 
test0 



For the Lboi of the two components, we use the results from lKing et all (|2010i ) of log(Lboi/£©) = 
—4.699 ±0.017 and —5.232 ±0.020. Their measurements come from integrating the resolved spectra 
of the two components from optical to mid-IR wavelengths 1^1 

The effective tempera tures of the two com ponents have been derived by fitting model atmo- 

(2009) analyzed resolved near-IR (JHK) spectroscopy of 
(|2006l ) model atmospheres computed with ±50 K spacing. 



spheres by two groups. (I ) iKasper et al 
the two components using 



Burrows et al 



They identified the best 3 fitting models to the data. We take the mean and RMS of their results, 
^eff A = 1275 ± 25 K and T e s b = 900 ± 25 K. The 25 K uncertainty is likely underestimated, as 
discussed in § 14.31 especial l y as t he fit to component Ba does not appear as good as that for com- 
ponent Bb. (2) iKing et al. ( 2010 ) obtained resolved 0.6-5.0 /iin spectra of the two components and 
compared the full spectra and the 1.25 /im K I absorption doublet to BT-Settl model atmospheres. 
Their comparisons applied an absolute flux calibration to the model atmospheres, based on the 



10 The two components of 2MASS J1209— 1004AB are subject to the exact same uncertainty in the photometric 
distance, and we take care to track this aspect properly when drawing the Monte Carlo values. The net result is that 
the uncertainty of the age difference between A and B is smaller than would be expected simply from the uncertainties 
in the absolute ages of either component. 

"Table |H indicates that 2MASS J22551861-5713056AB (q = 0.63), 2MASS J05185995-2828372AB (q = 0.54), 
2MASS J12255432-2739476AB (q = 0.63) have comparably small mass ratios as 2MASS J1209-1004AB and 
e Ind Bab. For these other binaries, the existing resolved measurements are insufficient to apply the isochrone 
test. 



12 



These values agree well with the Lboi derived from resolved Jf-band phot ometry and our bol ometric corrections 



in § 13.61 Using the high quality parallax to e Ind A (tt = 276.1 ± 0.3 mas; Ivan Leeuwenl 120071 ). the iGw/fO-band 
absolute magnitudes for the two components are Mk = 13.58 ± 0.02 and 15.84 ± 0.03 mag, and Tf-band bolometric 
corrections are 2.87 ± 0.10 and 2.20 ± 0.10 mag (including the intrinsic scatter and the ±0.5 subclasses uncertainty 
in the spectral types). Together, these give log^boi/^o) = —4.67 ± 0.04 and —5.31 ± 0.04 dex for e Ind Ba and Bb, 
respectively. 
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known luminosity and distance to the system. They found good agreement for T c g =1300-1340 K 
and 880-940 K for the two components. We adopt the mean values and use the quoted range as 
the RMS. 



Using the lKing et al.1 (120101 ) temperatures, Figure [6] shows the position of the two components 
of e Ind Bab on the H-R diagram, and Figure [7] shows the age difference between the two com- 
ponents, computed in the same fashion as § 14.11 Like for 2MASS J1209— 1004AB, the data are 
consistent with coevality for the binary, with a log(age) difference of 0.5^q' 3 and 0.3^04 dex from 
the Lyon and Tucson models, respectively, though only at the 2<r level when using the Lyon/COND 
models. The uncertainties are smaller for e Ind Bab co mpared to 2MASS J 1209— 1004AB due to 
the much more precise distance and temperatures. The iKasper et al.1 (|2009l ) temperatures give a 
similar result for the component ages, in fact slightly more consistent with coevality, though not 
significantly so given the uncertainties. 

Our conclusion of coevality for the system agrees with the (somewhat different) approach of 
Kasper et al.l (|2009l ). who determine log (g) and T e g- by fitting model atmospheres to the resolved 
spectra and t hen u se the results to compute the ages from evolutionary models. However, note that 
Smith et al.l (|2003l ) determine a hotter T e g of 1500 ± 100 K for e Ind Ba based on high resolution 
near-IR spectroscopy; if this temperature is adopted, then this component's H-R diagram position 
is off the model loci, with an age older than 10 Gyr and inconsistent with coevality. 

The age we derive from the H-R diagram is consistent with independent const raints on the age 



from t he primary star e Ind A. The most recently used estimate is 0.8-2.0 Gyr from lLachaume et al 



(|1999l ). based on their semi-empirical relation between stellar rotational period and age. Ap- 
pendix O updates the age estimate, leading to a broader range of 0.5-7.0 Gyr. Both values for 
e Ind A agree well with the median ages of 0.6-2.5 Gyr derived by us from the H-R diagram analysis. 



4.3. Systematic Uncertainties in T c fj's for T Dwarfs from Model Atmospheres 

For both 2MASS J1209— 1004AB and e Ind Bab, while the models do successfully indicate 
that the systems are coeval, the mass information derived from the H-R diagram reveals some 
interesting issues. Table [5] gives the computed component masses and mass ratios, based on the 
same calculations used to determine the ages. 

At face value, the HR diagram-derived mass ratio for 2MASS J 1209— 1004AB appears to be 
inverted, where the much brighter component A is non-sensibly indicated to be the lower-mass 
component. However, the formal uncertainties are very large and at the ~ la level, consistent 
with the q ~ 0.5 value inferred from evolutionary models (§ 13.71 and Table [3]). For e Ind Bab, the 
m ass ratio from the H-R diag ram is in good agreement with the value of q = 0.60 ± 0.02 derived 



by iMcCaughrean et al.1 (J2004J) , based on the component luminosities, evolutionary models, and an 
assumed age of 0.8-2.0 Gyr (where we have taken the average and RMS of their computed values). 
So within the accuracy afforded by the current data, the mass ratios derived from the H-R diagram 
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seem reasonable. 



However, a significant disagreement occurs when considering the total mass of e Ind Bab. The 
total mass derived from the H-R diagram is 72+g° Mj up (78^Jq Mj 1i n ) as derive d from the Lyon 



(Tucson) models and the iKasper et al.l (|2009l ) temperatures, with the King et al.1 (120101 ) tempera- 



tures giving a slightly higher total mass of 86±io M Ju P (95±i2 M Ju P )- (Note that the calculation of 
the total mass from the H-R diagram is independent of the aforementioned uncertainty in e Ind A's 
ag e). Such values strong ly conflict with the preliminary dynamical mass of 121 ± 1 Mj up measured 
by ICardoso et al.l (|2009l ) from the binary's relative orbital motion. And thus, there is a problem 
in one or more of the input data or models involved in the analysis. The observations of the 



compo nent Lboi's and total dynamical mass seem secure. The Lboi's are determined bv lKing et al 



(|2010l ) using resolved spectroscopy across a very broad wavelength range, agree well with the values 
obtained using near-IR photometry with our well-determined bolometric corrections (§ 14. 2p . and 
the uncertainties are fully propagated through our analysis. Likewise, the Cardoso et al. dynamical 
mass also seems reasonable given the large angular separation between the two components, the 
large number of observing epochs, and the long duration of the monitoring (~40% of the derived 



Bouv et al. 


— , . 
2004 


and 


Dupuv et al. 


2009a) 



Assuming for the moment that the loci of evolutionary models are correct, this discrepancy in 
total mass of e Ind Bab must arise from errors in the temperatures assigned to the binary compo- 
nents from the model atmosphere fittingly Generally speaking, changes to the relative tempera- 
tures between the two components would alter the inferred mass ratio, whereas systematic shifts to 
the temperatures assigned to both components would alter the derived total mass. Figures [8] and [9] 
summarize the changes to the HR-diagram-derived properties of the e Ind Bab system (component 
ages, age difference, mass ratio, and total mass) if the T e g 's for the two components are changed in 
these three ways: relative changes in the T e ff's, overall shifts in the T e g's, and changes to the T e g 
of only one component. These calculations place the following constraints on the T e ff estimates: 



• The relative T e g difference between e Ind Ba and Bb appears to be consistent with most of 
the constraints, namely coevality, an estimated age of ~1 Gyr, and a mass ratio of ~0.6 
(Figure [8]). However, the total mass of the system is under-predicted by the evolutionary 
models using these T e fj 's. 

• The disagreement between the total mass from the H-R diagram and the orbital monitoring 



13 A more exotic explanation would be that e Ind Bab is an unresolved triple system, with one of the components 
being an unresolved, nearly equal-mass binary. This would increase the total mass of the system derived from the 
H-R diagram, since correcting for such binarity would reduce the Lboi of the component and lead to a higher mass 
for a fixed T e ff. However, such a change would drive the H-R diagram positions towards non-coevality (assuming the 
assigned T e ff's are correct). More significantly, the resolved absolute IR magnitudes of the t wo known components o f 
e Ind Bab are consistent with other T dwarfs of comparable spectral types (e.g., Figure 16 of lBurgasser et al.ll2006dn . 
and thus correction for any unresolved binarity would make the associated component underluminous. 
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disappears if both components are ~50-100 K hotter (Figure [9]), which still leaves the system 
coeval. However, an older age of ~6 Gyr then results. 

• Changing the T e g for only one of the components cannot simultaneously fulfill the require- 
ments of coevality and agreeing with the dynamically measured total mass (Figure [TO]) . 



Altogether, this analysis points to a small, but non-zero, systematic error of ~50-100 K in 
the published T e g values from model atmosphere fitting. This seems plausible given the known 
inadequacies of the model atmospheres, especially in the L/T transition region (§ |T]), and thus 
invoking errors in the evolutionary models is not compelling. The implied temperatures for the 
two components would then be ~1350-1400 K and ~950-1000 K for components Ba and Bb, 
re spectively. For com ponent Ba, this hotter temperature is consistent with the 1500 ± 100 K found 
by lSmith et al.1 (120031 ) and noticeably hotter than previous fitting of early-T dwarfs (~ 1100-1200 K; 
§ 14. 1|) . For component Bb, the implied hotter temperat ure agrees with that fou nd expected for 
mid/late-T dwarfs based on model atmosphere fitting bv lBureasser et al.1 (12006a ). where a linear 



fit to their results gives 1030 K for the T6 subclass (Equation 7 of 



Liu et al 



2008T ) 



The systematic T e g errors that we infer from the isochrone method are comparable to the 
estimated uncertainties in spectral fitting of field T dwarfs (§ 14.11 and references therein), though 
e Ind Bab provides a more stringent test of the models given this binary's abundance of obser- 
vational constraints compared to ordinary field brown dwarfs 1^1 Also, note that the x-axes in 
Figures [8]- [TO] show the full possible range of T e g changes; larger changes would displace one or 
both binary components off the model loci on the H-R diagram, leading to ages that are either far 
too young (<10 Myr) or too old (>10 Gyr). 

Finally, we point out that given the measured -Lboi's and dynamical mass for e Ind Bab, Figured] 
shows that the evolutionary models would lead to an inferred age of about 5 Gyr. In the same 
fashion, King et al. (2010) have inferred an age of 3.7-4.3 Gyr from the evolutionary models (see 
also Cardoso et al. 2009). This is at the older end of the age range inferred for the K5V primary 
star e Ind A (Appendix IC]) . though the age-dating methods are largely derived for <1 Gyr solar- 
type stars and are not well-calibrated for such old ages or for such a late- type primary. In fact, it 
may be that the age of the e Ind ABab system is most accurately determined from brown dwarf 
evolutionary models, though it will be difficult to test this hypothesis independently. 



14 A somewhat similar, though less precise, consistency check of the atmospheric model fitting relative to the 
evolutionary models comes from the derived ages of field and benchmark T dwarfs, where the fitted T e g and log(g) 
values combined with evolutionary models lead to reasonable age determinations of a few to several Gyr for single 
field objects (e.g. iLeggett et aljEooj . 120071 : Isaumon et al"]|2007l : feurgasser et alll2006ah . 
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5. Conclusions 

We have identified the T3 dwarf 2MASS J 1209— 1004 as a tight substellar binary, with an 
extreme mass ratio (q ~ 0.5) compared to previously known field brown dwarf binaries. The near- 
IR photometry combined with modeling of the integrated-light spectrum indicates spectral types 
of T2.0±0.5 and T7.5±0.5 for the two components. This newly discovered binarity lessens the 
utility of this object as the primary near-IR spectral-type standard for T3 dwarfs. We suggest 
SDSS J1206+2813 replacement. 

The highly unequal mass ratio of the system allows us to test the joint accuracy of evolution- 
ary and atmospheric models, by examining if the ages inferred from the H-R diagram position of 
the two components are consistent with coevality. Using the photometric distance to the system 
and latest available T e g's from model atmosphere fitting of T dwarfs, we find that the models 
successfully indicate that the two components of 2MASS J1209— 1004AB are coeval, with a dif- 
ference of log(o<7e) = —0.8 ± 1.3 (— l.Ol^'g) dex based on the Lyon (Tucson) models. A similar 
analysis of the T1+T6 binary e Ind Bab also supports the accuracy of the models, with the system 
being consistent with coevality at the 1— 2a level, depending on the choice of model atmosphere 



fits and evolutionary models. (The temperatures for e Ind Bab from lKing et al.ll2010l in combina- 
tion with the Lyon/COND models produce the most non-coeval results, with an age difference of 
log(age) = 0.5±g;| dex.) 

While the models succeed in respect to the isochrones, for e Ind Bab the total mass derived 
from the H-R diagram (~80 Mj up ) is discrepant with the measured dynamical mass of 121 ±1 Mj up . 
This is most simply resolved by assuming a systematic increase in T e g for the two components by 
~50-100 K warmer, i.e., an increase to ~1375 K and 950 K for components Ba and Bb, respectively. 
Such an error is not unexpected, given the known inadequacies of the model atmospheres in the 
L/T transition region, the quality of the spectral fit, and the general uncertainties in the procedures 
used to fit such models to observed spectra. Such a T c g increase leads to an implied age of about 
6 Gyr for the system, somewhat older than previous estimates. 

While the ~50-100 K errors are comparable to those typically cited in the atmospheric model 
fitting, it is worth emphasizing that e Ind Bab indicates that such errors are present (i.e.. non-zero), 
as the masses inferred from the H-R diagram are strongly discrepant with the dynamical mass, a 
result that is independent of the uncertainty in the age of the system. However, assuming the 
evolutionary models are basically correct, T e g errors much larger than ~100 K are not warranted, 
as the inferred properties from the H-R diagram would then be inconsistent with the underlying 
constraints of coevality and a mass ratio <1.0. Also note that the very precise T e g results from 
the T5+T5.5 mass-benchmark binary 2MASS J1534-2952AB and the M4+T8.5 age-benchmark 
system Wolf 940AB point to ~100 K overestimates of T P ff when model atmo spheres are used to fit 



the near-IR spectrum of T dwarfs (|Liu et al.ll2008l : iBurningham et al.1 120091 ) , similar in amplitude 



though opposite in sign to what we find here from e Ind Bab. For Wolf 940B, iLeggett et al.l (|2010al ) 



show that luminosity-constrained model fits to the near-IR spectrum of this T8.5 dwarf are poor, 
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with the synthetic spectrum J and -ff-band peaks too high and the -ftT-band peak too low. Most 
likely this is due to the known incompleteness of the molecular opacity line lists in this region. 

Thus, substantial errors in evolutionary and/or atmospheric models for T dwarfs do not ap- 
pear to be required, though it is not possible to fully reconcile all the model predictions with the 
observations, espe cially for e Ind Bab w here the measurements are exceptionally precise. It is also 
worth noting that iDupuv et al.l (|2009bl ) have found that the luminosities predicted by evolution- 
ary models at ~0.8 Gyr may be systematically underestimated by a factor of 2-3, based on the 
luminosities and dynamical masses of the L4+L4 substellar binary HD 130948BC. A simple factor 
of 3 boost to the model-predicted luminosities seems implausible — the model tracks on the H-R 
diagrams (Figure [6]) would then give implausibly old ages (>>10 Gyr). Such a luminosity error 
could be made compatible with the H-R diagram positions if the current effective temperatures 
for the T dwarfs were also systematically overestimated by >200 K, but such a large discrepancy 
seems implausible. 

Even stronger constraints on the models from the isochrone test will require smaller uncer- 
tainties in Lboi and T e fj. For 2MASS J1209— 1004AB, numerical tests using the measurements 
from § 14.11 indicate that a ~5% distance and ~25 K uncertainty in T e g are needed to measure a 
2a discrepancy from coevality. In absence of a parallax, our analysis is consistent with coevality, 
but the results on absolute ages in Table [3] ca nnot be taken as a firm constraint. The parallax 
measurement is within current capabilities (e.g. IVrba et al.ll2004h and is being pursued by us using 

20091 ) . but the required T e s accuracy is 



the Canada-France-Hawaii Telescope (e.g. IDupuv &; Liu 



a challenge to atmospheric models. Both 2MASS J1209— 1004AB and e Ind Bab have plausibly 
short orbital periods for dynamical mass determinations, which will significantly strengthen the 
future tests of the models. This is especially true in the case of e Ind Bab, given the additional 
independent constraint from the age estimate of the primary star. 

The immediate prospects of finding similarly extreme systems to further test the models in this 
fashion are limited, given that hundreds of nearby ultracool dwarfs have already been imaged at 
high angu lar resolution bu t yieldi ng only 2MASS J1209— 1004AB and e Ind Bab. As one alterna- 
tive path, iBurgasser et al.1 (|2008al ) have pursued the possibility of finding previously unrecognized 
extreme mass-ratio ultracool binaries by modeling of integrated-light spectra, suitable for special 
cases where the primary and secondary spectral types are sufficiently discrepant (though the re- 
solved measurements of the individual components needed to place them on the H-R diagram may 
be challenging). In the n ear-future, the ady ent of powerful new wide-field surveys, s uch as the Pan 



STAR RS optical survey ([Kaiser et al.ll2002l ) and the WISE mid-IR survey satellite (jMainzer et al 



20091 ), will boost the prospects for expanding the census of field brown dwarfs and thereby enable 
discovery of new binaries. A larger sample of extreme systems will allow the isochrone test to be 
applied over a wider range of temperatures and luminosities, thus challenging the models over a 
large span of brown dwarf masses and ages. 
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A. Updated Near-IR Bolometric Corrections for Ultracool Dwarfs 



We improve t he polynomial fit for the i f -band bolometric correction as a function of near-IR 

by u pdating the spectral types for the T dwarfs in 



spectral type from IGolimowski et al.l (12004a 



their sampl e to types from 
HN Peg B dLiu et alihoOTl : 



Burgasser et a 



Luhman et al 



20061)1): adding data for the T dwarfs HD 3651B and 



2007bl : LLeggett et al.ll2008l ): and constructing fits for the 



J and .ff-band bolometric corrections. We use a 6th-order polynomial fit to accurately capture 
the behavior of the J and ii-band data and to avoid artificially steep changes in the polynomials 
at the earliest and latest spectral types. The results of the new fits are shown in Figure [11] and 
Table [6j Note that the use of the revised spectral types for the T dwarfs leads to 30-50% smaller 
RMS about the fit compared to the original Golimowski et al. study. 

We also plot the BC yalues for two of the latest type objects known, the T8.5 dw arf Wolf 940B 



( Burningham et al 



2007 



Smart et al 



20091; iLeggett et al.ll2010al l and the T9 dwarf ULAS 0034-0052 (jWarren et al 



2010i ). From the published bolometric luminosities, MKO photometry, and 



distances, we compute their bolometric corrections in JHK, though these are not included in the 
polynomial fit. The values are notably different from even the T7-T8 objects, likely due to the 
sharply increasing fracti on of Lhoi emitted at mid-infrared wavelengths in the very latest known 



objects (e.g. Figure 1 of lLeggett et alJl2010bl ) and perhaps impr ovements in the methods used to 
derive Lboi for late-T dwarfs using model atmospheres since the IGolimowski et al.1 (|2004al ) study, 
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which assumed a Rayleigh-Jeans approximation beyond 5 /im. 



B. Homogenous Compilation of Substellar Field Binaries 



We determined or compiled from the literature the spectral types, distances, and luminosities of 
all ultracool field binaries with integrated-light spectral types of L4 or later or those with ancilliary 
evidence for being substellar. The available literature compilations are based on a broad mix of 
input data and calculation methods. To improve upon this, we assembled a more homogenous 
compilation of binary properties, as described below. The final results are in Table HI 

Photometry: For the integrated-light photometry, we preferentially used measurements on the 
MKO photometric system where available and 2MASS measurements for the remaining objects. 

To determine the resolved photometry, we compiled the highest precision JHK flux ratio 
measurements available from the literature and supplemented them with measurements from public 
data archives (HST, VLT, and Gemini) and our own unpublished Keck LGS AO data. In cases 
where the integrated-light photometry and flux ratio s were not on the same ph otometric system, 
we used conversions between MKO and 2MASS from lStephens fc Leggettl (|2004l ). 



Infrared d ata were not ava i lable for 5 binaries originally discov ered by HST using the F814VF 



bandpass (e.g. iReid et al.ll200ll ; iGizis et al.ll2003l ; LBouv et al.ll2003l ). For these cases, we estimated 
the if-band flux ratio from the opt ical results. To do this, we assumed AF814W= AI, which is 
accurate to ~ 0.02 mag according to lBouy et al.1 (|2003l ). To estimate the secondary's spectral type, 
we derived the relation between spectral type and Mj, based on th e 11 objects that w ere used in 
the Liu et al. (2006) faint relations and have I-band photometry in lDahn et al.l (120021 ): 



Mj 



14.83 + (6.396 x 10" 1 ) x SpT - (1.542 x 10~ 2 ) x SpT 2 



(Bl) 



where SpT = means L0, etc, using optical spectral types for the L dwarfs. The fit is valid from 
L0 to T7.5, and the rms about the fit is 0.25 mag. Then using the adopted primary spectral type 
(described below) and the observe d AI, we d e rived the secondary spectral type. With spectral types 
for both components, we use the iLiu et al.l (|2006l ) SpT-absolute magnitude relations to estimate 
AK for the binary. We accounted for the additional uncertainty introduced by this process in our 
analysis, and the affected objects are marked in Table [Has being correspondingly uncertain. 

Resolved spectral types and distances: When available, we adopt the individual component 
spectral types previously determined from detailed analysis, e.g., spectral decomposition. 

For the remaining binaries, we assumed that the primary spectral types and uncertainties 
were identical to the integrated - light spectral types. We adopted optical spectral types deter- 
mined on the iKirkpatrick et al.l (|1999l ) classification scheme for L dwarfs and infrared types on 



the iBurgasser et al.1 (l2006bl ) classification scheme for T dwarfs, except in two cases. One was 
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DENIS-P J225210.7-173013AB, which iKendall et al.l ((200J) classified as L7.5 based on a com- 
parison of the H and if-band spectra to those of optically typed L dwarfs. The other was 
2MASS J09201 223+3517429AB, which has highly discrepant integrated-light spectral types of L6.5 
in the optical (jKirkpatrick et al.ll2000l ) and TOp in the infrared. We adopted the infrared type in 
this case. 

To derive the secondary sp ectral types, we used the average of the "faint" and "bright" SpT- 
absolute magnitude relations of iLiu et al.l (|2006l ) along with the known or estimated absolute mag- 
nitudes for the secon daries. (The Liu et al. relations were extended to objects of M6-L1 using 
photometry from the lLeggett et al.ll2010bl literature compilation.) For binaries with parallax mea- 
surements, we computed the absolute magnitudes of the secondary components directly. For the 
rest, we estimated phot ometric di s tance s based on the primary components' apparent magnitudes 
and the same averaged ILiu et all (|2006l ) SpT-absolute magnitude relations. If distance estimates 
were available from multiple bandpasses, we preferred the if-band estimate, then i^-band, and 
finally J-band, to minimize the potential non-monotonic behavior in the L/T transition region. 
The uncertainties in the resulting distances account for the input photometry uncertainties and the 
intrinsic scatter in the empirical relations. 

Using the final assembled IR absolute magnit udes for the s econda ries, we estimated spectral 
types using the method described in Section 3.2 of iDupuy et al.1 ()2009cl ). For a given bandpass, we 
computed the spectral type probability distribution for an object based on its absolute magnitude, 
the uncertainties in both the photometry and the distances, and the intrinsic scatter in the empirical 
SpT-absolute magnitude relations. We then combined the results from all available bandpasses 
to form the joint spectral type probability distribution and thus to estimate the final secondary 
spectral types and uncertainties. We rounded all spectral types and their uncertainties to the 
nearest 0.5 subclass. 



Bolometric luminosities and luminosity ratios: We derived individual luminosities for the 
binary components from each available bandpass (JHK) using the resolved photometry, distances, 
and bolometric corrections. Bolometric corrections were computed based on the polynomial fits 
in Appendix lAl and the spectral types assigned to each component. For computing luminosities, 
we assume Mbd,© = 4.75 mag. Uncertainties for the resulting Lboi's were computed from the 
uncertainties in the resolved photometry, in the bolometric corrections due to the spectral type 
errors, and the intrinsic scatter in the bolometric correction relations. 

We also computed luminosity ratios for each available bandpass; this is typically known to 
higher precision than Lt, i since the ratios are independent of the distance uncertainty. The error 
in the luminosity ratio comes from the quadrature sum of the errors in the flux ratios and the 
bolometric corrections of the primaries and secondaries. For our final adopted luminosities and 
luminosity ratios, we chose the bandpass that gave the smallest error in the luminosity ratio. 



Mass ratios: We used the 



Burrows et al 



(119971 ) evolutionary models to derive the mass ratios for 
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each system, using the derived Lboi's, their uncertainties, and an adopted age. We take care to treat 
the distance uncertainty as an error common to each binary, as opposed to applying it independently 
to each individual component. For objects with specific age information in the literature, we use 
that; otherwise we assume 1 Gyr. In both cases, we ignore the age uncertainty in calculating 
the mass ratio, since this is common to (most of) the objects. Four binaries have their total 
mass measured dynamic ally, and we us e their published mass ratio measure ments in these cases: 
2MASS .T15 34-2952AB JLiu et alJ boOfih. HD 130948BC ([Dupuv et alJl2009bh . LP 349-25AB, and 



GJ 569Bab (|Dupuv et al 



201 



To summarize the mass ratio results, we fit the observed values with a power-law distribution 1^1 
Since the sample is relatively sparse, instead of the usual process of fitting to a binned histogram, 
we adopt a maximum likelihood approach. The probability distribution of the mass ratio q is 
described by 

P(q) = {a + 1) q a (B2) 

for g = — l, with the constant chosen so that the total probability is unity. We define the log of 
the likelihood function as 



C = lnP(g[a) 



\n\[(a + l)q? 
1=1 

n 

y~] [ In (a + 1) + a In 

=i 

n 

ln(a + 1) + a In q t 



i=l 



(B3) 
(B4) 
(B5) 



i=i 



where n is the number of objects in the sample and q% is the mass ratio for each object. To find 
the best fitting a, we find the maximum of C: 



n 



dC 



da a + 1 ^— ' 

i=i 



+ In qi 







(B6) 



which gives the final result 



a 



1 + 



n 



(B7) 



For our sample, the distribution C(a) is well- fit by a gaussian, from which we extract la confidence 
limits. To incorporate the uncertainties in q^, we use a Monte Carlo approach, computing the best 



15 We exclude the L6p+T7.5p system, SDSS J141624.08+134826.7AB, from this ensemble analysis given its unusu- 
ally large separation and mass ratio compared to all other field systems. Also, because of its large Lboi ratio, its mass 
ratio estimate is highly dependent on the assumed age, from 0.2 to 0.5 for ages of 1 to 10 Gyr. 
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a many times with the qi values drawn from normal distributions; the resulting RMS dispersion in 
a from the ensemble of results is much smaller than the 68% confidence interval in any realization 
of L. As a simple approximation, we add the RMS (cr a = 0.3) in quadrature to the confidence 
limits (cr a = 0.6) to obtain our final uncertainty on a. 

We find a = 4.9 ± 0.7 for the sample in Table [H For such an exponent, binaries with q < 0.65 
such as e Ind Bab and 2MASS J1209— 1004AB amount to only 9% of t he population. Our resu lt is 
comparable to, though somewhat steeper, than the g( 4 - 2±L0 ) found bv lBurgasser et al.l ( 2006dl ) for 
a sample of 30 binaries with estimated primary masses of <0.075 Mq. There are some differences 
in the fitting method and the input data; for the latter, the Burgasser et al. compilation included 
young binaries in open clusters and star-forming regions (which tend to have a broader range of 
mass ratios) and used the heterogenous set of mass ratios published in the literature. Burgasser et 
al. also corrected the observed mass ratio distribution to account for the bias towards equal-mass 
systems in mag nitude-limited surveys , based on an approximate conversion between flux ratio and 
mass ratio (see iBurgasser et al.ll2003l ). We choose not to apply any corrections, since the surveys 
involved in Table S] span a wide range of wavelengths and telescopes. Given these differences, a 
more shallow power-law is expected from the Burgasser et al. comp i lation . Similarly, using a full 
Bayesian analysis to account for differing survey sensitivities, lAllenl (|20071 ) derived a = 1.8 ± 0.6; 
his sample included many late-M (stellar) primaries which likely led to a shallower exponent. 



C. Summary of Age Estimates for e Ind A 



Age determinations for the K5V star e Ind A are necessarily indir e ct an d rely on empirical 
correlations between stellar activity, rotation, and age. lLachaume et al.l (|1999l ) estimated a mean 
log(age) = 9.1 (1.3 Gyr) with a range of 8.9 - 9.3 (0.8 - 2.0 Gyr), based on their semi-empirical 
relati on between st e llar r otational period and age and the estimated rotational period of e Ind A 
from lSaar Sz OsterJ ( 1997 ) as derived from Ca II measurements and an empirical relation between 
chromospheric acti yity and rotation. This age range has been used for analysis o f e Ind Bab in 
previous work (e.g. iKasper et al.ll2009i ; iMcCaughrean et al.l 12004 ; IScholz et al.ll2003l ). 



Additional age estimates come from the correlation bet ween age and chro mospheric activity 
traced by Ca II line emission, as measured by the R'uk index. iHenry et al.l ()1996l ) measured e Ind A 
on three occasions spanning 1.3 yr, with a mean and an rms of logi?^^- = —4.56 ± 0.01 dex. 



Chromospheric activity is known to be variable on decade-long time scales, and lGray et al.l ((2006) 
found log R' HK = —4.85 for e Ind roughly eight years later. 



The available age calibrations of R'uk give consiste nt estimates. ( 1) Us ing thejDonahuc 



calibration, we find ages o f 1.2 Gyr and 3.2 Gyr fro m the lHenry et al.1 (j 19961 ) and lGray et al 



1998) 



20061 ) 



data, respectively. (2) The lSoderblom et al.l (1 199 if) calibration results in ages of 1.2 Gyr and 3.4 Gyr. 
(3) Note that with a B — V color of 1.06 mag dPerrvman et al.lll997h. e Ind A is beyond the most 
recent activity-age calibrations, those of iMamajek <fe Hillenbrand! (120081 ) which are only valid for 
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0.5 < B — V < 0.9 mag stars. If we disregard this limitation, their R' HK calibration gives ages of 
0.9 and 4.0 Gyr for the two set s of measurements. (4) We can also use the "preferred" method of 
Mamaiek fc Hillenbrand! (|2008l ). namely converting R' HK to a Rossby number, in order to estimate 
a rotation period that can then be used to derive an age from their gyrochronology relation. The 
resulting ages are 1.6 Gyr and 5.0 Gyr for the higher and lower activity levels, respectively. The 
nominal uncertainty for such chromospheric ages is esti mated to be around ±50% at a single epoch 



(Soderblom et al 



1991 



Mamajek &: Hillenbrand! 120081 ) . Thus, the overall age range derived from 



chromospheric indicators is 0.5-7.0 Gyr, with a mean of about 2 Gyr. 



An alternative estimate comes from lJanson et al.l (|2009l ). They obtained high-contrast imaging 
of the system to attempt to directly image a close-in substellar companion, w hose existence ha s 
been inferred from a long-term linear trend in the radial velocity of e Ind A (jEndl et al.1 l2002l ) . 
By combining their imaging non-detection combined with the radial velocity data, Janson et al. 
conclude that the star is likely to be older than than ~1 Gyr, or else the companion would have 
been sufficiently self-luminous be have been directly detected. 

This age estimate may be improved in the future by measuring e I nd A's rotati o n p eriod, as 
gyrochronology relations can pr ovide more precise age estimates (15-25%; 



Barnes 



rotati o n p eriod, as 
sl l2007h R Based on 



Mamajek &; Hillenbrand (j2008l ) , the rotation period of e Ind A is expected to be 14-50 days given its 
level of chromospheric activity. Furthermore, because e Ind A is a very bright star (V = 4.7 mag), 
it is within reach of asteroseismological measurements with current facilities, which may be able to 
provide an even more precise age estimate. 
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Fig. 1— Keck LGS AO imaging of 2MASS J1209-1004AB. North is up and east is left. Each 
image is 0.75" on a side. Contours are drawn from 90%, 45%, 22.5%, 11.2%, and 5.6% of the peak 
value in each bandpass. 



- 32 - 



O 



X 
I 



1.0 - 



0.5 - 



0.0 - 



-0.5 - 



-1.0 - 




-1.0 



-0.5 0.0 
(H-K) MK0 



0.5 



Fig. 2. 
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Warren et a 



Birmingham et al 



I2004) and 


Chiu et al. 


Pinfield et al. 


f 


2008) 



( 2008J) : and in d ividual compone n ts of re so lved binaries fromlMcCaughrean et al 



Liu et al 



|200J), iBureasser et alJ (|2005bl 1. iBurgasser et all (j2006dh . iLiu k Leggettl (|2005l ). and 
(|2006l ). (For the T8.5 dwarf CFBDS 0059—0114, the if -band measure ment is synthesized from 



the near-IR spectrum and i^s-band photometry of iDelorme et al.ll2008l .) The photometry errors 
are comparable to or smaller than the size of the plotting symbols for most of the objects. The 
numbers indicate the near-IR spectral subclass of the objects, with half subclasses being rounded 
down (e.g., T3.5 is labeled as "3") , and objects of the same subclass plotted in the same color. 
The late-L dwarfs (classified on thelGeballe et al.ll2002l scheme) are plotted as bare numbers. The 
T dwarfs (on the IBurgasser et al.1 l2006bl scheme with extension by iBurningham et al .1 120081 ) are 
plotted as circumscribed numbers, with squares for integer subclasses (e.g., T3) and diamonds for 
half subclasses (e.g., T3.5). 
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Fig. 3.— The integ rated-light near-IR spectrum of 2MASS J1209-1004AB (heavy black line; 



Burgasser et al.l l2004l ) modeled as the sum of two T dwarfs (green lines): one early-T template 
(blue lines) and one late-T template (red lines). The templates are chosen based on their simi- 
lar JHK colors to the two components of 2MASS J 1209— 1004AB. The agreement between the 
observed spectrum and the modeled blends is good, meaning that the spectral types inferred for 
2MASS J1209— 1004A and B from the near-IR photometry are consistent with the integrated- light 
spectrum. 
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Fig. 4.— Inferred T cS and log(#) for the comp onents of 2MASS J12 09-1004AB and e Ind Bab 
based on solar-metallicity evolutionary models by iBurrows et al.l (|l997l ) . The blue-colored hatched 
regions show the constraints from the observed Lbol's of the binary components and the inferred 
system ages, with the A component for each system being the region on the right (i.e., the hotter 
object). For 2MASS J1209-1004AB, we plot a generic age range of 0.1-10 Gyr. For e Ind Bab, we 
plot an age range of 0.5-7.0 Gyr, with the filled circle representing 2 Gyr. The uncertain -Lboi's of 
2MASS J1209— 1004AB are reflected in its relatively wide hatched regions; these can be improved 
with a parallax determination for the system. The solid lines are isochrones from 0.1 to 10 Gyr, 
labeled by the logarithm of their age; the two unlabeled isochrones are 0.05 and 0.1 Gyr. Dotted 
lines represent iso-mass models, labeled on the right of the plot in units of Mq. The vertical 
grey regions show the T e g's for the two co mponents derived fro m model atmosphere fitting; for 
e Ind Bab, we sho w the T e g values from both lKasper et al.l (|2009l ) (cooler values, in light grey) and 
King et al.l (|2010l ) (hotter values, in dark grey). 
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Fig. 5. — Mass ratios for all known field brown dwarf binaries, based on our compilation given 
in Tableland described in Appendix [Bl The mass ratio for 2MASS J1209— 1004AB from our 
calculations is shown as the filled color bar, demonstrating the rarity of such an unequal mass ratio 
among the field population. (SDSS J141624. 08+134826. 7AB is not plotted, given its very unusual 
properties relative to all other known systems and the large systematic uncertainty in its mass 
ratio; see Appendix IBl) 
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Fig. 6. — H-R diagram showing the two components of 2MASS J1209— 1004AB and e Ind Bab 
compared to isochrones and iso-mass tracks from the Tucson and Lyon (COND) evolutionary 
models. The data points show the individual components. The thick black diagonal connects the 
two components, showing the isochrone defined by the binary system. The isochrones are shown as 
solid colored lines, with the numbers in squares on the left side of the plots giving the logarithms 
of the isochrone ages (10 7 ,10 8 ,10 9 ,10 10 yr). The iso-mass tracks for 0.01 M & (upper line) and 
0.04 Mq (lower line) are shown as dashed black lines. For 2MASS J1209-1004AB, two sets of 
errors are shown for the y-axis: the longer, thinner grey error bars show the total error in the L^ \ 
measurements, and the shorter, thicker black errors bars show the errors with the uncertainty in 
the distance removed, since t his uncerta i nty is common to both components. For e Ind Bab, the 
plot uses the T e g values from iKing et al.l (|2010l ) . 
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Fig. 7. — Probability distribution of the age difference between components A and B in the 
2MASS J1209— 1004AB and e Ind Bab systems, based on ages inferred from their locations on 
the H-R Diagram (Figure [6]) and the Lyon (COND) and Tucson evolutionary models. The his- 
togram shows the results from our Monte Carlo calculation, and the thin solid line shows the best 
fitting Gaussian to the binned probability distribution. The dashed vertical line denotes coevality 
of the two components. 
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Fig. 8. — The effect of changing the relative T e g values assigned to e Ind Ba and Bb on the 
physical properties derived from the H-R diagram positions, based on the Lyon/COND evolutionary 
models. The relative T e ff difference between the two components is assumed to be applied with an 
equal but opposite amplitude to both components. For instance, the leftmost point represents a 
50 K decrease in T e g for component Ba and a 50 K increase for component Bb. The four plotted 
properties are the ages of the individual components, the age difference, the mass ratio, and the 
total mass. The horizonal dashed lines represent the available constraints: for this system, the age 
of the primary star (log(i/yr) = 9.1), the requirement of coevality (log^/ts) = 1), the mass ratio 
inferred from the luminosity ratio (q = 0.6), and the total d ynamical mass (121 Mj up ) . The grey 
points are the results based on using the Kasper et al. ( 20091') TLff's a nd th e green on the King et al. 
(|2010l ) T e g's. The rightmost point is missing for the iKing et al.l (|2010l ) models as such a large 
temperature change moves the binary components off the loci of available models. 
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Fig. 9. — The effect of ap plying a const a nt off set to the assum ed T P ^ value s for e Ind A and e Ind B, 
with grey being r e g from lKasper et al.l ()2009l ) and green from lKing et al.l (|2010l ). For example, the 
leftmost point represents a 150 K decrease from the nominal T e g values employed in § 14.11 See 
Figure [8] caption for more details. 
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Fig. 10. — The effect of changing the assumed T e g- value for one co mponent of e Ind Ba b while 
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and green from iKing et al.l (|2010l ). Left: Changing the T e g- of component Ba. For example, the 
leftmost point represents a 100 K decrease from the nominal T e g value employed in § 14.11 See 
Figure [8] caption for more details. Right: Changing the T e g of component Bb. 
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Fig. 11. — Bolometric correction as a function of near - IR spe c tral type for th e MKO J, H, and 



if -band niters, based on d a ta from [ Golimo wski et a 



(|20n7h . lLeggett et all (j^nnsl ). ISmart et al.1 ( 20ld ). and lLeggett et al.1 (j2ni0ah . Points without visible 



(12004ah. iLuhman et alJ (|2007bh . iLin et al 



error bars have uncertainties smaller than the plotting symbol. The solid line overplots a 6th-order 
polynomial fit to the datasets for objects as cool as T8, with c oefficie nts given in TableEl The dotted 
line for the if-band data gives the fit from lGolimowski et al.l (|2004al ). Known and suspected binaries 
are plotted as open circles and as circumscri bed small c ircles , respectively, and excluded from the 
fits. Suspected binaries are those flagged by iLiu et al.1 (|2006l ) based on their large over lum i nosity 
relative to objects of similar spec tral type. The L dwarfs are classified on the lGeballe et al.l (J2002) 
scheme and the T dwarfs on the iBurgasser et al.l (|2006bl ) scheme. The two coolest objects (T8.5 
Wolf 940B and T9 ULAS 0034-0052) are excluded from the fit. 



Table 1. Keck LGS AO Observations 



Date 


Filtcr a 


Airmass 


FWHM 


Strchl ratio 


Separation 


Position angle 


Amag 


(UT) 






(mas) 




(mas) 


(deg) 




2007-Apr-22 


./ 


1.46 


112 ± 3 


0.022 ± 0.002 


143 ± 6 


311 ± 2 


1.5 ± 0.2 




H 


1.40 


105 ± 5 


0.049 ± 0.004 


151 ± 13 


314 ± 5 


2.8 ± 0.3 




K s 


1.53 


92 ± 5 


0.121 ± 0.019 


156 ± 7 


314 ± 4 


3.2 ± 0.5 


2008-Jan-16 


CH 4 s 


1.22 


77 ± 5 


0.083 ± 0.012 


142 ± 7 


313 ± 3 


2.07 ± 0.17 



Note. — All photometry on the MKO system. The tabulated uncertainties on the imaging performance 
(FWHM and Strehl ratio) are the RMS of the measurements from the individual images. 



Table 2. Resolved Properties of 2MASS J1209-1004AB 



Property Component A Component B 



J (^magj 


15.80 ± 0.05 


17.27 ± 


1 7 


CH4S (mag) 


15.20 ± 0.04 


17.27 ± 


0.15 


H (mag) 


15.32 ± 0.04 


18.08 ± 


0.26 


K s (mag) 


15.12 ± 0.14 


18.29 ± 


0.53 


K (mag) 


15.22 ± 0.04 


18.48 ± 


0.52 


J — H (mag) 


0.48 ± 0.06 


-0.81 ± 


0.31 


J — K (mag) 


0.58 ± 0.07 


-1.21 ± 


0.54 


CH 4 s - H (mag) 


-0.12 ± 0.05 


-0.81 ± 


0.30 


H - K (mag) 


0.10 ± 0.05 


-0.40 ± 


0.58 


Estimated spectral type 


T2.0 ± 0.5 


T7.5 ± 


0.5 


dphot (pc) 


21 


± 4 




M(J) (mag) 


14.19 ± 0.41 


15.66 ± 


0.44 


M(CH 4 s) (mag) 


13.59 ± 0.41 


15.66 ± 


0.44 


M(H) (mag) 


13.71 ± 0.41 


16.47 ± 


0.49 


M(K S ) (mag) 


13.51 ± 0.43 


16.68 ± 


0.67 


M(K) (mag) 


13.61 ± 0.41 


16.87 ± 


0.66 


\og(L boLA+B /L Q ) 


-4.61±0.15(0.02) a 




Mboi (from J-band) 


16.21 ± 0.45 


18.24 ± 


0.47 


log(L bo i/L G ) (from J-band) 


-4.58 ± 0.18 


-5.39 ± 


0.19 


\og(L boLA /Lboi,B) (from J-band) 


0.81 ± 0.11 





Note. — All infrared photometry on the MKO photometric system. 
We also measure intcgrated-light Spitzer/IRAC photometry for the sys- 
tem of [3.6] = 14.02 ± 0.03, [4.5] = 13.49 ± 0.03, [5.8] = 13.33 ± 0.03, 
and [8.0] = 13.06 ± 0.07 mags (§ EU). 

a The error in parenthesis gives the Jboi error assuming there is no 
uncertainty in the distance. 
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Table 3. Properties of 2MASS J 1209- 1004 AB Derived from Evolutionary Models 



Property Component A Component B 



t = 0.5 Gyr 



Total mass (M Jup ) 52 ± 7 

Mass (M Jup ) 36 ± 6 17 ± 4 

q(=M B /M A ) 0.47 ±0.16 

T e g(K) 1281 ± 162 771 ± 103 

T cSyB /T cff , A 0.60 ± 0.08 

Radius (R Q ) 0.101 ± 0.002 0.109 ± 0.003 

log(g) (cgs) 4.96 ± 0.11 4.56 ± 0.13 
Orbital period (yr) 27+^° 



t = 1 Gyr 



Total mass (M Jup ) 75 ± 10 

Mass (M Jup ) 49 ± 9 25 ± 5 

q (= Mb/M a ) 0.51 ± 0.10 

T eff (K) 1336 ± 173 802 ± 107 

0.60 ± 0.07 

Radius (R Q ) 0.092 ± 0.002 0.101 ± 0.003 

log(ff) (cgs) 5.18 ± 0.11 4.80 ± 0.12 

Orbital period (yr) 22±| 2 



t = 5 Gyr 



Total mass (M Jup ) 132 ± 10 

Mass (M Jup ) 77 ± 5 55 ± 8 

q (= M b /Ma) 0.72 ± 0.07 

T eff (K) 1418 ± 163 881 ± 120 

T cSB /T cSA 0.62 ± 0.07 

Radius (i? Q ) 0.082 ± 0.001 0.083 ± 0.003 

log(g) (cgs) 5.48 ± 0.04 5.31 ± 0.11 



Orbital period (yr) 17 +31 



t = 10 Gyr 



Total mass (M Jup ) 149 ± 8 

Mass (Mjup) 79 ± 3 70 ± 7 
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Table 3 — Continued 



Property Component A Component B 

q (= Mb /Ma) 

Teff{K) 
Teg,B/T e S,A 

Radius (Rq) 
log(g) (cgs) 
Orbital period (yr) 



0.90 ± 0.03 
1424 ± 159 908 ± 119 
0.64 ± 0.07 
0.081 ± 0.002 0.078 ± 0.002 
5.50 ± 0.03 5.47 ± 0.09 
15.9tf 8 4 



Table 4. Homogenous Compilation of Substellar Field Binaries 



Binary 


SpT(A+B) SpT(A) 


SpT(B) 


d (pc) 


Ago (Gyr) 


iogf-tbol.A) 


lo g(£bol,B) 


, , L hal,A , 
1 ° S( I^ ) 


9 
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005 


40 


18 


26 


61 


18 


a 


Gl 417BC 


L4 


5 


± 





5 


L4 


5 


± (1 


5 


T , - + 10 

L4.5„ 5 


21.9 ± 0.2 


0.74 (0.58-0.86) 


-4 


06 ± 0.03 





03) 


-4 


18 + 


(14 





04) 





12 


= 





(15 





94+0 


02 


3 


16, 42 


82 


61 


2 




GJ 1001BC 


L5 





± 





5 


LB 





± 11 


5 


L5.0 ± 1.0 


13 0+°' 7 
1JU -0.6 


3 (1-10) 


-4 


03 ± 0.06 





051 


-4 


08 + 


07 





05) 





05 


± 





08 





99+0 


01 


25 


16, 25 


33 


20 


38 


e 


DENIS-P J1228.2-1547AB 


LB 





± 





5 


L5 





± 11 


5 


L5.0 ± 1.0 


20 - 2 -o'.7 




-4 


19 ± 0.05 





031 


-4 


24 + 


or, 





04) 





05 


± 





05 


I) 


98+0 


03 


35, 45 


16, 42 


31 


17 






2MASS J12392727+5515371AB 


LB 





± 





5 


L5 





± (1 


5 


L6.0 + ?S 
— 1.5 


19 ± 4 




-4 


23 ± 0.17 





04) 


-4 


31+0 


17 





06) 





00 


± 





07 


(1 


95+0 


05 


24 


3, 16 


32 






b 


2MASS J21321145+1341584AB 


L6 


o 


+ 





5 


L5 





± (1 


5 


L7.5 ± 0.5 


23 ± 4 




-4 


39 ± 0.17 





04) 


-4 


72 + 


17 





04) 





33 


± 





04 


(1 


76+0 


08 


58 


16, 58 


14, 58 








2MASS J17281150+3948593AB 


L7 









5 


L7 





± 


5 


T Q n+ 20 
L9.0 _j 


24 ± 2 




-4 


47 ± 0.08 





04) 


-4 


65 + 


08 





04) 





18 


± 





04 


(1 


85+0 


08 


24 


16, 42 


32 


62 






2MASS J08503593+1057156AB 


L6 





i 





5 


L(i 





± 


5 


L9.5 i 1.0 


26± 3 2 




-4 


53 ± 0.08 





03) 


-4 


87 + 


0!) 





03) 





34 


± 





04 





75+0 


03 


50 


37, 42 


31 


17 






2MASS J22551861-5713056AB 


LG 


(J 


-4- 


1 





LB 





± 1 





T,8 fl + 1 n 


12 ± 3 




-4 


38 ± 0.20 





04) 


-4 


95 + 


22 





07) 





57 


± 





00 





63+0 


05 


53 


9, 16, 51 


53 








2MASS J21522609+0937575AB 


L6 


o 


± 


1 





L6 





± 1 





L8 5+ 4 '° 
° -3.0 


22 ± 5 




-4 


36 ± 0.19 





04) 


-4 


39 + 


20 





04) 





03 


± 





05 





98+0 


04 


51 


16, 51 


53 








SDSS J141624.08+134826.7AB 


L6 


o 


-j- 





5 


L6 





± (1 


5 


T7.5 ± 0.5 


8±l 




-4 


30 ± 0.20 





03) 


-6 


24 + 


20 





03) 


1 


94 


± 





04 





18+0 


01 


12 


12 


5, 12 


57 




f 


DENIS-P J225210.7-173013AB 


L7 


5 


± 


1 





L4 


5 


± 1 





T4.5 ± 1.0 


11 ± 2 




-4 


51 ± 0.18 





04) 


-4 


93 + 


18 





04) 


(1 


42 


± 





05 





70+0 


03 


52 


9, 16, 42, 51 


11, 30 






g 


2MASS J05185995-2828372AB 


Tl 





± 





5 


L7 


5 


± 11 


5 


T5.0 ± 0.5 


20 ± 4 




-4 


55 ± 0.19 





10) 


-5 


29 + 


23 





15) 





75 


± 





22 


I) 


54+0 


10 


9 


9, 16 


10, 11 








SDSS J042348. 57-041403. 5AB 


TO 





± 


(1 


5 


L7 


5 


± (1 


5 


T2.0 ± 0.5 


15.2 ± 0.4 




-4 


33 ± 0.04 





03) 


-4 


59 + 


04 





03) 





26 


± 





04 


(1 


80+0 


03 


7 


9, 37 


10, 11 


62 






DENIS-P J020529.0-115925AB 


L7 





± 





5 


L7 





± (1 


5 


L6.0 ± 1.0 


19.8 ± 0.6 




—4 


28 ± 0.04 





04) 


—4 


34 + 


05 





04) 





06 


± 





05 





96+0 


03 


35 


16, 42 


31 


17 






2MASS J09153413+0422045AB 


L6 





± 


1 





L(i 





± 1 





L7.0±f° 


18 ± 5 




-4 


28 ± 0.20 





04) 


-4 


33 + 


21) 





04) 





05 


± 





05 


I) 


97+0 


04 


51 


9, 16, 42, 51 


53 








2MASS J21011544+1756586AB 


L7 


5 


± 





5 


L7 


5 


± 11 


5 


T1.0± 2 2 ° 


33+4 




-4 


55 ±0.11 





04) 


-4 


67 + 


13 





08 





13 


± 





09 


I) 


89+0 


07 


3 


3, 13 


32 


62 




b 


2MASS J03105986+1648155AB 


L8 





± 





5 


L8 





± (1 


5 


L9 -°ll'.6 


25 ± 4 




-4 


54 ± 0.15 





03) 


-4 


60 + 


15 





04) 





06 


± 





04 


(1 


95+0 


08 


50 


37, 42 


32 








Gl 337CD 


TO 





± 





5 


TO 





± 


5 


T0.5ti 5 


20.4 ± 0.2 


1.5 (0.6-3.4) 


-4 


62 ± 0.04 





04) 


-4 


70 + 


or, 





05) 





08 


± 





04 





93+0 


03 


8 


16, 42 


10, 63 


61 


63 




2MASS J09201223+3517429AB 


TO 





± 





5 


TO 





± 


5 


T1.0±%-° 


18 ± 3 




-4 


75 ± 0.15 





04) 


-4 


84 + 


15 





04) 





00 


± 





05 





93+0 


04 


50 


16, 42 


10, 32 






h 


SDSS J105213. 51+442255. 7AB 


TO 


5 


± 


1 





TO 


5 


± 1 







22 ± 4 




-4 


76 ± 0.16 





03) 


-4 


78 + 


10 





04) 


(1 


112 


± 





05 





98+0 


04 


42 


13, 42 


13 








SDSS J205235. 31-160929. 8AB 


Tl 





± 


1 





Tl 





± 1 





Tl K + l ■ ^ 

io.o_ 2 


25 ± 5 




-4 


74 ± 0.16 





03) 


-4 


85 + 


10 





04) 





12 


± 





05 





90+0 


04 


42 


13, 42 


13 








2MASS J14044948-3159330AB 


T2 


5 


± 


(1 


5 


TO 


(1 


± 


5 


T5.0 ± 0.5 


18 ± 3 




-4 


84 ± 0.16 





05) 


-4 


99 + 


10 





05) 





16 


± 





06 





88+0 


04 


44 


44 


11, 43 








SDSS J102109. 69-030420. 1AB 


T3 





± 





5 


Tl 


I) 


± 


5 


T5.5 ± 0.5 


29+\ 




-4 


58 ± 0.13 





04) 


-4 


75 ± 


13 





05) 





17 


± 





06 





86+0 


or, 


9 


9, 37, 42 


10, 11 


60 






e Ind Bab 


T2 


5 


± 


11 


5 


Tl 


11 


± 


5 


T6.0 ± 0.5 


3.622 ± 0.004 


2 (0.5-7.0) 


-4 


62 ± 0.03 





03) 


-5 


24 + 


04 





04) 





(il 


± 





05 





63+0 


02 


48 


16, 48 


10, 48, 56 


61 


1 




SDSS J153417.05+161546.1AB 


T3 


5 


± 





5 


Tl 


5 


± 


5 


T5.0 ± 0.5 


31 ± 6 




-4 


87 ± 0.16 





03) 


-5 


09 + 


10 





04) 





23 


± 





05 


I) 


83+0 


03 


40 


13, 40 


40 








2MASS J12095613-1004008AB 


T3 





± 


11 


5 


T2 


11 


± 


5 


T7.5 ± 0.5 


21 ± 4 




-4 


58 ± 0.18 





07) 


-5 


39 + 


10 





00) 





81 


± 





11 





51 + 


05 


1 


1, 13 


1, 10 








SDSS J092615. 38+584720. 9AB 


T4 


5 


± 


(1 


5 


T4 


5 


± 


5 


T3 '°-4.0 


22 ± 4 




-4 


79 ± 0.16 





06) 


-4 


95 + 


21 





08) 





Hi 


± 





13 


(1 


88+0 


00 


9 


9, 37 


10 








2MASS J15344984-2952274AB 


T5 





± 


(1 


5 


T5 





± 


5 


T5.5 ± 0.5 


13.6 ± 0.2 


n _„+0.09 


-5 


02 ± 0.02 





02) 


-5 


09 + 


02 





02) 


(1 


08 


± 





02 


0.936+0 


012 


6 


34, 41 


10, 41 


60 


41 


a 


2MASS J12255432-2739476AB 


T6 





± 





5 


TO 


I) 


± 


5 


T7.5 ± 0.5 






-4 


93 ± 0.05 





05) 


-5 


47 + 


or, 





04) 





54 


± 





06 


0.63+0 


03 


6 


37, 42 


10 


60 






2MASS J15530228+1532369AB 


T7 





± 





5 


T7 





± 


5 


T7.5 ± 0.5 


11 ± 3 




-5 


42 ± 0.24 





04) 


-5 


54 + 


25 





03) 





12 


± 





05 


0.90+0 


05 


9 


9, 34, 42 


10 









Note. — Field binaries arc included in this table if the primary spectral type is L4 or later, or if there is other evidence that the system is young enough to be substellar, i.e., from accretion signatures 
(DENIS-P J0041-56AB), lithium absorption (Kelu-1 and 2 MASS J1146+22AB), low surface gravity (DENIS-P J0357-44AB and SDSS J2249+00AB), or dynamical masses (LP 349-25AB and Gl 569Bab). Distance 
references arc listed for objects with trigonometric parallaxes or those with a wide main-sequence companion with a good distance estimate; otherwise the distances arc photometric estimates derived by us, except 
DENIS-P J0041 — 56AB comes from Rcincrs (2009). The uncertainties on £<bol are listed with and without the uncertainty in the distance folded in. Other notes: (a) The mass ratios come from the published 
analysis using the measured total dynamical mass and bolomctric luminosity of the system combined with the luminosity ratio of the two components, (b) The .RT-band flux ratio was estimated from _F814W-band 
flux ratio (Appendix B). (c) The age estimate assigned by us is based on the presence of lithium and the £<bol of the components using evolutionary models (e.g., Liu & Lcggctt 2005). (d) This object also has 
a distance estimate from a comoving wide F8 companion star (G 171-58) of 42.2^^'g pc (Faherty et al. 2010), which is somewhat discrepant with our photometric distance listed here, (e) GJ 1001BC has an 
integrated light spectral type of L5.0i0.5 and essentially identical fluxes and colors; thus both components arc likely to be L5 ( Golimowski ct al. 2004b). The absolute magnitudes of component C formally given 



to a spectral type of L4.0ztl.O using our methodology in Appendix^] which is consistent. We adopt L5.0rfc0.5 for both components, (f) Burningham ct al. (2010) and Burgasscr ct al. (2010b) find the secondary is 
high gravity and low temperature and suggest an age of ^10 Gyr and 2—10 Gyr, respectively. For an age of 10 Gyr, the mass ratio of the binary would be ftiO.5, based on their mass estimates for component B and 
that of[Bowlcr ct al. ( 201(^) for component A. (g) We list the IR spectral type, as an optical type is not available for this object, (h) This has very disparate integrated- light spectral types, L6.5 in the optical and 
TO in the near-IR. We use the IR type here. 
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Table 5. Ages and Masses Inferred from the H-R Diagram 



1 eff,A^ 1 eff,B IVLOQei 


iog(Wy r ) 


log(Wyr) 


\og(t A /t B ) 


M A 
(M Jup ) 


M B 
(M Jup ) 


Mb + Ma 
(M Jup ) 


9 = M s /M^ 


2MASS J1209 1004AB 


1150 ± 100 K, 800 ± 70 K (this work) Lyon/COND 

Tucson 


7 5+ 1 - 2 

' -°-1.5 
' -°-1.5 


8.8j# 


-i-oiij 


ioi 7 e 

11+24 
ii -8 


ist% 

20±t? 


o 4 +27 
^-19 
o fi +63 


i y+5.2 

1 6+ 6 ' 2 


e Ind Bab 


1275 ± 25 K, 900 ± 25 K (Kasper et al. 2009) Lyon/COND 

Tucson 

1320 ± 20 K, 910 ± 30 K (King et al. 2010) Lyon/COND 

Tucson 


9.1±0.2 
9.1±0.2 
9.4±0.2 
9.3±0.2 


8.8±0.3 
9.0±0.2 
8.9±0.3 
9.0±0.3 


0.2±0.3 
0.1±0.3 

n 5 +0 - 4 

u -°-0.3 

3 +0 - 3 

U.O_ 4 


49±8 
57±7 
63l? 


27±£ 
28±? 
30±8 


72+i° 

86±i 
95t» 


0.6±0.2 
0.6±0.2 
52+ - 17 
n 4Q+ 019 



Note. — The results for each component are the median and 68% confidence limits from the Monte Carlo-generated probability distributions 
for the ages and masses computed from the H-R diagram positions (§HJ). 



Table 6. Coefficients of Polynomial Fits for Bolometric Correction Versus Near-IR Spectral Type 



Bandpass 


CO 


Cl 


C2 




C3 




C4 




C5 




C6 




RMS (mag) 


Jmko 


1.890448c+01 


-8.053993C+00 


1.491738e- 


fOO 


-1.367605e- 


01 


6.540717e- 


03 


-1.558986e- 


04 


1.462266e 


-06 


0.14 


Hmko 


1.366709e+01 


-5.426683c+00 


1.027185e- 


fOO 


-9.618535e- 


02 


4.733874e- 


03 


-1.171595e- 


04 


1.148133c 


-06 


0.07 


Kmko 


5.795845c+00 


-1.471358e+00 


2.868014c- 


-01 


-2.647474e- 


02 


1.282657e- 


03 


-3.177629e- 


05 


3.159780e 


-07 


0.08 



Note. — These are the coefficients of the 6th-order polynomial fits plotted in Figure [TT] for the bolometric correction as a function of 
near-IR spectral type. The fit for filter X is defined as 

6 

BC{X) = Y, Cl x (SpTY 

i=0 



where the numerical spectral type is defined as SpT =11 for LI, Sp T = 12 for L2, SpT = 20 for TO, etc. The fits are applica ble from 
M6 to T8.5. The L dwarfs are classified on the iGeballe et al.l (|2002[ ) scheme and the T dwarfs on the iBurgasser et al.l (]2006b| ) scheme. 
The last column gives the standard deviation about the fit in magnitudes. 



